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Abstract 

Starting  from  the  differential  equation  of  the  continuity  of  the  current 
flow^  a  general  solution  is  given  for  probleiis  in  which  the  conductivity  is 
given  as  a  function  of  space  and  time.  A  physical  Interpretation  for  the 
conQ;>lete  Maxwell  current  is  obtained^  and  it  is  shown  how  the  Maxwell  current 
can  be  eongiosed  of  two  field  vectors^  namely,  that  of  the  electrostatic  field 
and  of  the  stationary  current  flow.  Each  vector  is  weighted  by  a  different 
time  function,  which  can  be  calculated  from  the  time  function  of  the  current 
source. 

A  method  is  developed  for  the  calculation  of  the  stationary  current  flow 
field;  and  the  eigen  functions  for  cartesian,  cylindrical,  and  polar  coordi¬ 
nates  are  given.  The  mirror  lav  for  the  ciurrent  flow  Is  determined  in  a 
medium  where  the  conductivity  increases  with  altitude  according  to  an 
e-funotion.  The  mathematical  formalism  is  explained,  using  the  exao^le  of  the 
field  of  a  decaying  point  source. 
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ON  a!HE  TBEOBY  OF  THE  ATMOSPHERIC  ELECTRIC  CURRENT  FLON,  17 


INTRODUCTION 

Problems  in  atmospheric  electricity  are  usTially  treated  according  to  the ' 
electrostatic  theory  in  spite  of  the  fact  that  a  three-dimensional  current 
flow  is  being  dealt  vith  here.  The  reason  for  this  unsatisfactory  state  is 
that  the  electrostatic  theory  is  worked  out  in  great  detail  in  numerous  text 
books  of  loathematicB  and  physics^^^^'^  but  problems  of  the  three-dimensional 
current-flow  theory — especially  with  a  variable  conductivity  or  a  convection 
cimrent— are  only  briefly  discussed,  if  at  all.  Therefore,  the  author 
started  to  publish  some  of  his  notes  on  the  current-flow  theory  under  the 
general  title,  "Zur  Stroemungs-theorie  des  luftelektrischen  Feldes  I,  II, 

III."  Also,  two  papers, which  deal  with  the  calculation  of  the 

three-dimensional  current  flow  under  certain  bovindary  conditions. 

Preceding  the  author,®  Holzer  and  Saxon^  published  their  famous  paper, 
"Distribution  of  Electrical  Conduction  Currents  in, the  Vicinity  of  Thunder¬ 
storms,"  and  in  1953  Taimira  presented  at  the  First  International  Conference 
on  Atmospheric  Electricity  his  excellent  "Analysis  of  Electric  Field  after 
Lightning  Discharges."  In  both  papers  the  calculation  was  carried  out 
according  to  the  current-flow  theory.  But  in  the  majority  of  publications  in 
the  field  of  atmospheric  electricity,  the  electrostatic  theory  is  still  used 
for  calculations,  and  this  quite  often  leads  to  erroneous  conclusions. 

It  soon  became  obvious  that  the  solution  of  a  few  single  problems  was 
not  enough  to  introduce  the  application  of  the  theory  of  current  flow  to 
atmospheric  electric  problems  in  general.  A  more  systematic  discussion  of 
the  methods  of  solution,  with  a  strong  accent  on  the  physical  meaning  of  the 
mathematical  manipulations  and  the  consequences  of  the  introduced  assun^tlons, 
is  necessary.  This  discussion  is  given  in  this  report. 

DISCUSSION 

The  Differential  Equation  of  the  Three-Dimensional  Current  Flow 

Table  I  gives  a  list  of  the  symbols  used,  with  their  meanings  and 
dimensions. 

In  the  electrostatic  theory,  only  the  quantities  on  the  left  side  of  the 
list  are  present;  while  in  the  current-flow  theory  there  are,  in  addition, 
the  fovir  quantities  of  the  rig^t  side.  It  Is  seen  that  the  dimensions  of  the 
quantities  of  a  corresponding  pair  differ  only  by  the  addition  of  the  time 
xinit,  8,  to  the  dimension  of  the  electrostatic  quantity.  This  will  prove 
very  fortunate  I'n  the  conversion  of  solutions  of  the  electrostatic  theory 
into  these  of  the  corresponding  solution  of  the  current-flow  theory,  as  will 
be  shown  later.  However,  this  "fortunate"  relationship  has  a  deeper  physical 
meaning,  which  becomes  much  clearer  if  the  dimension  of  the  charge  cb 
(coulomb)  Instead  of  the  current  A  (as^re)  is  introduced  into  the  dimension 
of  the  quantities. 


Table  I.  Symbols  Usedi  vlth  Their  Meanings  and  Diinenalons 


current  density 


Q  [ab]  ■  charge  I  [a]  «  strength  of  current  source 

q,j^  •  space  charge  density  ”  space  current  source  density 

^  ~  «  dj.8leetrie  displacement  _r  A-i 

L  nr J  1  ~  B  current  density 

r-Ae-i 

e  r,.  dielectric  constant  tat 

J  "  strength 

^  [v]  tt  potential  function 


With  the  relation  As  »  cb,  the  foUoving  are  obtained: 


ft  [eb] 

^  m 

^  [$\ 

“  [^] 

”  [f] 

'  [a 

r  cbl 

[vmj 

^  [is] 

In  this  way  the  time  has  been  removed  from  the  electrostatic  quantities, 
where  it  does  not  belong,  and  introduced  into  the  current-flow  theory,  where 
one  would  e:^ct  it  to  be.  If  the  time  is  now  extended,  e.g.,  one  second  to 
infinity,  then  the  movement  of  all  particles  would  freeze  and  everything 
would  be  static.  The  fo\ir  specific  quantities  of  the  current-flow  theory, 

I,  (u,  i,  and  would  vanish  because  the  time  is  in  the  denominator  of  the 
dimensions,  and  the  electrostatic  conditions  are  obtained.  So,  one  could  say 
that  electrostatic  is  a  special  case  of  the  current-flow  theory. 

In  the  current-flow  theory,  all  the  quantities  listed  above  can  be  func¬ 
tions  of  space  and  time.  The  problems  encountered  are  mostly  of  the  follow¬ 
ing  type:  Given  are  the  spatial  distribution  and  the  time  function  of  the 
current  source  is  and  the  electrical  properties  of  space  e  and  \  as  functiozis 
of  space  and  time.  This  Includes  also  the  boundary  condition  such  as  a  per¬ 
fect  conducting  earth  surface,  or  sometimes  an  ionosphere  of  infinite  con¬ 
ductivity.  The  problem  is  to  calculate  the  potential  function  and  the  field 
and  current  distribution. 
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The  dielectric  constant  e  is  alvays  a  true  constant  in  space  and  timej 
since  the  difference  of  e  in  air  and  in  a  vacuum  is  negligible.  The  manner 
in  vhich  the  conductivity  X  is  given  determines  the  grade  of  diff ic\>lty  of 
the  problem^  and  to  some  e^rtent  also  the  method  of  solution.  Therefore,  the 
problems  may  be  subdivided  into  the  foUoving  six  classes: 

1.  The  conductivity  is  constant  in  space  and  time. 

2.  The  conductivity  is  constant  in  space,  but  a  function  of  time. 

3.  The  conductlylty  is  constant  in  time,  but  a  function  of  space. 

4.  The  conductivity  is  a  function  of  space  and  time. 

The  conductivity  is  a  tensor  (Influence  of  the  earth's  magnetic 
field  in  the  ionosphere). 

6.  The  assuz^tion  for  the  definition  of  the  conductivity  is  no  longer 
valid. 


This  is  the  case  in  outer  space,  vhere  the  mean  free  path  of  the  moving 
particle  is  greater  than  either  the  antenna  of  the  measuring  instrument  or 
the  space  under  consideration. 

Classes  1,  2,  and  3  are  special  eases  of  class  4.  Class  4,  again,  could 
be  considered  as  a  special  case  of  However,  the  solution  of  ^  as  veil  as 
of'  6  class  problems  requires  a  such  larger  soathematlcal  effort.  Therefore, 
this  report  is  limited  to  a  general  solution  of  1  to  4  class  problems. 


The  solution  to  class  1  problems  can  be  obtained  easily  if  the  corres* 
ponding  electrostatic  solution  is  knovn.  It  is  necessary  only  to  substitute 

for  the  electrostatic  quantities,  Q,  q,  S,  and  e,  the  corresponding  current- 
flow  quantities  I,  (u,  iT,  and  X  of  table  I  in  the  fomulas  of  the  potential 
function  0  or  of  the  electric  field  distribution  £.  For  instance,  the 
potential  function  ^  i  of  a  point  charge  leads  to  the  potential  func- 

T  1  ^ 

tlon  0  ■  ~  ^  point  current  source. 


The  applicability  of  class  1  problems  in  atmospheric  electricity  seems 
to  be  very  limited  as,  for  Ixuitance,  to  the  calculation  of  the  effective 
altitude  or  area  of  field  or  current  antenzmis,  where  the  conductivity  in  the 
space  considered  may  be  assumed  to  be  constant.  But  it  is  shovn  later  that 
the  electrostatic  solution  is  needed  in  addition  to  that  of  the  stationary 
current  flow  for  the  general  solution  of  class  3  and  4  problems.  Also,  the 
iz^rtant  matching  condition^  for  air-earth-current  mesisurements  can  be 
obtained  with  the  assuz^ion  that  the  conductivity  is  constant  in  space  and 
time. 


If  the  conductivity  changes  with  tizie  but  is  constant  in  space,  this 
would  result  in  mismatching,  which  is  a  problsz  belonging  to  class  2.  An 
analysis  of  these  conditions  for  alr-earth-current  measurement  at  the  ground 
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l8  given  by  Ruhnke.  At  a  larger  scale  tbe  same  problem  is  encountered  by 
the  air-earth-current  radiosonde, ^3  where  the  conductivity  increases  daring 
the  flight  time  from  the  low  ground  value  by  a  factor  100  or  more  to  that  in 
the  higher  altitudes.  Here  the  eissujnption  is  made  that  the  conductivity  can 
be  considered  as  constant  in  the  space  occupied  by  the  antennas  of  the  radio¬ 
sonde,  but  vlU  Increase  In  time  vlth  the  ascent  of  the  balloon. 

!Ihe  calculations  given  by  the  author^»'^»^^  belong  to  class  3» 

There  is  no  exangile  of  class  4  problems  in  the  literature.  This  type 
of  problem  will  be  encountered  by  the  calculation  of  the  austausch  generator, 
where  the  air  turbulence  vlU  change  the  conductivity  in  space  and  time  and 
also  cause  a  convection  current  which  is  a  function  of  space  and  time. 

The  differential  equation  of  tbe  three-dimensional  current  flow,  which 
has  to  be  fulfilled  by  all  solutions ,  is  derived  from  the  second  Maxwell 
equation  and  is  well  known  as  the  equation  of  the  continuous  current  flow. 

It  is  written  in  the  following  fonn: 


div 


at 


+  i) 


w. 


(1) 


The  intensity  of  the  ctirrent  production  is  usually  given  by  a  convection 
current .  Q.  In  this  case  the  space  current  source  density  u)  is  given  by 

div  Q.  IS  this  is  introduced  into  Eq.  (l), 


div  (^  +  i)  »  div 


From  Eq.  (2)  follows 


SD  7  -r  -* 

^  +  1  »  ;  +  C. 


(2) 

(3) 


c,  which  appears  here  as  a  kind  of  integration  constant  with  the  condition 

div  c  =  0,  is  known  as  the  con^ilete  Maxwell  current.  It  is  a  current  flow, 
which  originates  at  the  boundaries  of  the  considered  space.  In  thunderstorm 

problems,  c  would  represent  the  fair-weather  current  between  the  ionosphere 
and  the  earth,  which  penetrates  the  tinonderstora  areas  as  well  as  the  fair- 
weather  areas.  However,  the  fair-weather  current  density  is  much  smaller 
than  the  current  density  caused  by  the  thunderstorm,  and  therefore  it  is 
usually  neglected.  This  cannot  be  done  in  the  case  of  the  austausch  genera¬ 
tor,  which  worlss  in  the  fair-weather  region.  Here,  the  convection  cnxrrent 

density  of  the  austausch  generator  Q  and  the  current  c  produced  by  the 
charged  earth  are  of  the  same  order  of  magnitude.  (The  conduction  current 

density  of  is  the  well-known  oir-earth-current  density.)  Therefore,  in 
this  case,  both  current  sources  have  to  be  considered.  But  as  the  principle 
of  superposition  also  holds  in  the  current-flow  theory,  the  current  dis¬ 
tribution  of  the  two  sources  may  be  calculated  separately,  and  the  solutions 
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superin^osed.  The  dielectric  displacement  D  and  the  current  density  i  are 
connected  vlth  the  electric  field  E  through  the  following  two  equations: 

e  ?  »  eT. 

X  E  «  1 .  (5) 

(If  X  Is  a  tensor  Aas  iQ  class  5  problems,  A  *  F •*  T,  and  from  here  on  the 
calctilation  of  class  5  problems  would  branch  off.)  Vith  Eq.8.  (4)  and  (5), 

E  is  substituted  for  T  and  D  in  (3),  and 

(6) 

is  obtained.  The  solution  of  this  differential  equation  is 


f 


,^(t^)  +  ^  J  (T  +  c) 


(7) 


This  is  the  general  solution  for  W  for  all  class  1  to  4  problems.  The  inte¬ 
gration  constant  gives  the  field  distribution  for  t=0.  On  account  of 

the  integration  constant,  the  integral  can  be  written  without  boxmdaries;  and 
the  letter  t,  which  ordinarily  would  appear  as  the  upper  boundary  of  the 
integral,  is  used  as  the  integration  variable  of  the  integreO.. 


According  to  the  above-mentioned  principle  of  superposition,  the  field- 
vector  E*’may  be  split  up  into  two  parts;  namely,  which  results  from  the 

convection  current  X>  and  ©b,  which  results  from  ciunrent  emitted  at  the 
boundaries . 


(8) 


and 


Eb 


e 


(9) 


Equation  (8)  cannot  be  discussed  without  further  Isnowledge  of  the  con¬ 
vection  current  "g.  The  analysis  of  Eq.  (9)^  especially  the  physical  meaning 

of  the  conqplete  Maxwell  current  ?  is  given  below. 


The  Solution  of  Boundary  Problems  without  Convection  Currents 

Set  ^  »  0,  and  from  Eq.  (6)  the  electric  field  Eb  due  to  the  boundary 

conditions  may  be  obtained.  With  the  further  omission  of  the  subscript  b, 
it  is 
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®  +iis 

at  e 


(10) 


The  form  and  the  physical  meaning  of  the  complete  Maxwell  current  c  vill 
now  he  deduced  from  the  three  boundary  values  c^,  c^,  and  ^  which  "c  has  to 
assume,  first  at  the  boundary  (^),  second,  if  X  =  0  (Fg),  and  third,  if 

f  -»  t's)- 

At  the  boundary  of  the  current  source,  the  field  Et,  is  given  by  the 
field  strength  and  the  time  function  T  of  the  current  source.  Therefore, 


Eb 


(11) 


F'j,^x,y,s)  1®  hereby  a  function  of  the  space  coordinates  x,  y,  and  z 

only,  and  is  a  pure  time  function.  This  in^ilies  that  on  all  parts  of 

the  boundary  of  the  current  source  the  field  changes  with  the  same  time 
function.  From  Eqs.  (lo)  and  (ll)  follows 


e 


^h  dt 


T. 


(12) 


If  X  =  0,  the  electrostatic  field  distribution  is  obtained,  which  will  be 
indicated  by  a  subscript  e  at  the  Maxwell  current  "oq  and  the  field  vector 

P^.  Furthermore,  tlie  finite  propagation  velocity  of  the  electric  field  with 
the  speed  of  ligiit  will  be  neglected,  which  means  that  the  electric  field 

E  follows  everywhere  tVje  time  function  T  of  the  source  momentarily.  E  = 

Fg  •  T.  Hereby  excltulod  are  all  problems  where  the  propagation  velocity  of 
electric  signals  becomes  impor-OTt,  for  instance,  by  the  electromagnetic 
Wave  of  a  lightning  dj}3oha,rge  (sferics).  From  Eq.  (lO)  follows 


(13) 


If  the  time  function  of  the  current  source  Is  a  constant,  the  condition  of 
the  stationary  current  flow  is  obtained,  which  is  indicated  by  the  subscript 

s.  With  =  0  and  ¥3  Fg  T,  from  Eq.  (10)  is  obtained 

^  ^  T.  (14) 

e  e  ^ 


Equations  (12),  (I3),  and  (l4)  indicate  the  form  which  the  l»laxwell  cur¬ 
rent  ?  in  general  will  assume.  It  may  be  inferred  that  "c  is  composed  in  the 
following  manner: 


e  F, 


dt 


+  X  T. 


(15) 
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Equation  (15)  fulfilJ.s  all  the  required  conditions.  At  the  botindary, 
vhere  Fgi^  «=  Fgi^  =  F^,  Eq.  (I5)  changes  to  Eq.  (12),  and  therefore  meets  the 
boundary  conditions.  Furthermore,  dlv  =  0,  because  div  Fg  =  0  as  well  as 
div  =  0.  Tlierefore  the  complete  l-laxwe].!  current  c  is  given  in  the  form 

of  Eq.  (15). 


The  j)hyGlcal  meaning  of  Eq.  (I6)  beco.iies  much  cleai'er  if  the  following 
fCbreviations  are  introduced  for  the  time- function  factors  of  Fg,  and  Fg. 


_  p  Xdt 
'■  e 

f 

Tg  =  e 

J  dt 

®  dt. 

(17) 

.  r 

J  e 

‘  Xdt 

Tg  =  e 

•'  e 

®  dt. 

(18) 

_  r 

To  =  e 

e 

• 

(19) 

By  partial  integration  of  Eq.  (17)> 

T  =T-e  e 

e  J  e 

is  obtained,  and  with  Eq.  (I8), 

Te  =  T  -  Tg,  T  =  Tg  +  Ts.  (20) 

The  sum  of  the  time  functions  of  the  electrostatic  field  and  the 
stationary  current-flovr  field  is  the  time  function  of  the  current  source. 
These  time-function  factors  are  obtained  by  the  prescribed  integration 
process  on  the  time  function  T  of  the  source.  They  are  different  at  each 
point  in  space,  because  X  is  a  function  of  spacej  and  they  are  of  different 

form  for  the  electrostatic  field  Fg,  the  stationary  current-flow  field  Fg, 
and  the  decaying  initial  field 

If  Eqs.  (17),  (18),  and  (I9)  are  introduced  into  (16) , 
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(21) 


r «  f;  +  %  Tg  +  f;  t^. 

The  electric  field  E  is  composed  of  three  vectors,  each  of  them  modulated  hy 
a  different  time  function.  Fq  is  the  field  distribution  at  the  time  t  »  0 

and  decays  with  the  time  constant  of  the  observation  point.  is  the  elec¬ 
trostatic  field,  and  the  wei^tlng  function  Tg  is  zero  for  a  constant- 
current  source,  but  becomes  very  large  for  fast  time  variations  of  the 

source.  is  the  field  of.  the  stationary  current  flow,  and  the  factor  Tg 
becomes  one  for  a  constant  current  source.  It  is  pointed  out  that  in  general 

and  Fg  do  not  have  the  same  directions,  and  the  rules  of  vector 
addition  have  to  be  applied.  After  the  decay  of  T^,  the  field- vector  F  is 
limited  to  the  space  between  Fg  and  Fg. 

To  illustrate  the  point,  the  field  lines  of  a  point  soiree  are  drawn  in 
Fig.  1  for  the  electrostatic  ease  (broken  lines),  and  for  the  stationary 
current-flow  case  in  which  the  conductivity  increases  with  altitude  according 
to  an  e-functiou  (solid  lines).  It  is  seen  that  for  almost  any  point  in 
space  the  field  direction  and  amplitude — the  amplitude  is  indicated  by  the 
spacing  of  the  lines- -Is  different  in  the  two  cases.  If  one  arbitrary  obser¬ 
vation  point  is  selected,  a  picture  of  the  composition  of  the  field-vector 
E  is  obtained,  as  given  in  Fig.  2.  The  wel^ting  time  functions  will  change 

with  time,  and  accordingly  the  vector  F  will  change  Its  amplitude  and  shift 
its  position  inside  the  hatched  space. 

For  a  fast  time  variation,  F  is  predominantly  given  by  F ;  and  for  slow 

Hr  ^ 

time  veurlatlon,  E  will  move  back  to  the  position  of  Fg.  Hence  it  becomes 
obvious  that  for  a  complete  solution  the  electrostatic  field  has  to  be  cal¬ 
culated  as  well  as  the  field  of  the  stationary  current  flow.  The  electro¬ 
static,  the  current  flow,  and  the  initial  field  Fg,  end  Fq  are  given  by 

the  gradient  operation  from  their  potential  functions  and 

respectively. 

fe  -  -  «rad  0*,  •  -  grad  0»,  F^,  -  -  grad  0,,.  (22) 

Hence,  the  electric  field  E  is  alto  given  by  the  gradient  of  a  potential 
function  0.  F  »  -  grad  0.  Therefore,  with  respect  to  Sq,.  (21), 

grad  0  *  To  grad  0©  ♦  T*  0i-  (23) 

Equations  (21)  and  (23)  ere  the  general  solutions  to  all  class  1  to  4 
problems  in  the  absence  of  a  eorvection  current. 

Solution  of  the  mfferential  Squatjon  of  the  Stationary  Current  Flow 

If  the  current  flow  is  stationary,  derivatives  with  respect  to  time  are 
zero,  and  from  Eq.  (2)  is  obtained  the  differentiel  equation  of  the  station¬ 
ary  current  flow. 
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Fig.  2 


.  Conrposition  of  the  generail  field  vector  E  "by  vector  addition 
from  the  electrostatic  field  vector  vei^ted  by  the  time 
function  Tg  and  the  stationary  current  flow  field  vector  ?g 
veighted  by  the  time  function  Tg. 
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(24) 


div  ig  =  0. 

With  ifi  B  X  5^  and  fg  «=  -  grad  0g, 

div  X  grad  0s  “  0*  (25) 

The  following  self-explanatory  operations  are  carried  out  on  this  equation; 
div  X  grad  0  *=  0. 


X  div  grad  0  +  grad  X  grad  0*0. 

Divide  hy  \^,  add  and  subtract  0  div  grad  X^. 

1.  1.  ^ 

X2  div  giad  0  +  2  grad  X^  grad  0+0  div  grad  X^  -  0  div  grad 

xi  «  0. 


div  (xl  0)  -  0  div  grad  (X^ 


0. 

0. 


(26) 


Here  are  defined  two  new  functions; 

M  a  X2  0  and  H  =  X^.  (27) 

Introducing  these  new  fvinctions  into  Eq.  (24),  the  relatively  simple 
equation  is  obtained; 


With  this  equation  there  are  available  aQJL  the  solutions  of  the  electro¬ 
static  theory  div  grad  M  =  0,  if  it  is  possible  to  present  X  by  a  function 

of  which  div  grad  N  *  0. 

It  is  pointed  out  that  different  coordinate  systems  may  be  used  for  the 
representation  of  M  and  N  because  the  operation  div  grad  is  inveuriant 
against  a  change  of  coordinates. 

As  an  exanQple,  assume  that  it  is  desired  to  ceuLculate  the  potential 
function  of  a  point  source  Imbedded  in  a  mediian  where  the  conductivity  can 
be  represented  or  approximated  in  the  space  xmder  consideration  by  a  suitable 
piece  of  a  parabola. 

X  »  Xq  (az  -  a)^.  (29) 

The  parameters  m,  a,  and  Xq  can  be  used  to  select  the  best-fitting  piece  of 
the  parabola. 

From  Eq.  (27)  follows 

i  1  , 

H  =  X^  =  X^5  (mz  -  a). 
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(30) 


and 


div  grad  N 


(31) 


div  grad  N  «=  =  0. 

dz^ 

Introducing  this  result  in  Eq.  (28)^ 

div  grad  M  ■=  0.  (32) 

This  equation  means  that  any  known  electrostatic  potential  function  or  any 
sum  of  them  may  be  chosen  In  any  coordinate  system  which  suits  the  purpose 
best.  For  a  point  source,  polar  coordinates  would  certainly  be  chosen, 
where  the  function  M  is  ^ven  by 

M  =  A  .  (33) 


A  is  here  a  constant  given  by  the  strength  of  the  source  I,  and  r  is  the 
distance  of  the  obseirvatlon  point  from  the  center  of  the  point  source.  If 
the  point  source  is  located  In  the  cartesian  coordinates  at  the  point  Xq, 

yo,  Zq,  then  r  is  given  by  r  -  [(x*Xq)^  +  poten¬ 

tial  function  is  now  easily  obtained  by  Eqs .  (27)»  (29),  and  (33) • 

0  -  ^  - •  (34) 

(m?.  •  a)  r 

The  determination  of  A  remains  as  the  last  step.  If  the  point  source  is 
inclosed  by  a  small  sphere  of  the  radius  Tq,  which  is  chosen  so  small  that 
the  conductivity  X  nan  be  considered  as  constant  with  the  value  at  the  cen- 

p 

ter  of  the  sphere  X  «=  A<}  (s>  •  a)  ,  then,  according  to  Eq.  (34)> 


Xq^  (m  -  a)  r 


The  field  in  the  r  direction  on  the  little  sphere  r  *  is  given  by 


Xo^  (m  -  a)  Tj, 


and  tlie  current  density  according  to  01m*B  law  i  *  XP, 

2  2 

A  Xq  (m  So  -  a)  A  Xq  (»  Zq  " 

i  a  . . . . . . . —  a  - 

^1  2  2 
Xq^  (m  Zo  -  a)  Tq 


Integration  over  the  surface  0  of  the  little  sphere  gives  the  current 
strength  I  of  the  source. 


I  =  J  ir  dO  = 


(^ 


4  ^ 
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or 


I 


(36) 


4  TT  (b  Zq  -  a) 

If  this  result  is  introduced  in  Eq..  (34),  the  final  solution  for  the  poten¬ 
tial  function  of  the  problea  is  obtained. 


im  Zrt  -  aj 


(D  z  -  a)  •  r 


Unfortunately,  in  atmospheric  electric  problems,  the  conductivity  is 
usually  represented  by  an  e-function.  For  instance,  in  a  cartesian  coordi¬ 
nate  system  x,  y,  z,  in  the  manner 

X  -  Xo  (38) 

The  x,y  plane  with  z  «  0  represents  hereby  the  earth's  surface,  and  the 
conductivity  at  the  ground.  This  function  (given  In  Eq.  (38))  does  not 
fulfill  Laplace's  equation,  but  It  is 


dlv  grad  V 


m  k®  e^  «  k^'N, 


dlv  grad  H_ 


Therefore,  the  function  M  is  given  as  the  solution  of  the  differential 
equation 

div  grad  M  -  k^  M  -  0.  (l»0) 

To  solve  Eq.  (40),  the  method  of  the  eigen  functions  viU.  be  applied, 
but  confined  to  problems  of  rotational  symmetry.  This  means  that,  for 
instance,  in  polar  coordinates  r,  9,  0,  the  function  M  is  given  by  the 
product  of  a  function  f(r);  vhieh  depends  on  r  only,  and  another  function 

g^Q^,  vhieh  depends  oxily  on  9. 

M  ■  ^(r)  ®f0)"  (^X) 

This  leads  with  Eq.  (4o)  to  the  differentiad  eqtiation 
dlv  grad  f ’g  -  k^  f ‘g  «  0, 


or,  in  polar  coordinates. 


d^  ^  2  df 


+  ^  Bln  0  ^  -  k2  fg  .  0, 

Bln  9  ^ 


d%  ^  2  df 
r  dr 


r^  k2  g  sin  9  as 


i-  sinO  ^-1*0.  (43) 


or 


1 

k^  f 


If 


g  sin  e 


ae 


sin  9 


ae 


(n+l) 


m 


where  n^O,  1,  2,  3^  •  •  •>  the  spherical  functions  of  the  first  p 
(Legenare’s  functions)  ea^  obtained  as  a  solution  for  g,  ^ 


®  “  ^n(cos  O)”  (^5) 

If  Eq.  (44)  is  inserted  in  (43)  ,,  the  differential  equation 

5  ^  /  B.  f  =  0  is  obtained  for  f .  (46) 

^  ^  \  r2  / 


Equation  46  can  be  solved  with  the  help  of  cylinder  functions  with  a 
noninteger  parameter  p  --  n  +  The  solution  is  given  by  ^ 

f  .  (j  k  r)"2  (J  k  r).  (47) 

Thereby,  j  is  the  imaginary  imit  ^ 


The  con^ilete  solution  for  M  follows  from  Eqs.  (4l),  45),  and  (47). 

Mn  *  An  Pn  (5  4  r)"?  (j  k  r) .  (48) 

Ajj  is  an  arbibr8.i7r  constajit,  which  serves  to  meet  the  prescribed 
boundary  eondltio.nu  M  may  be  presented  in  a  more  general  way  as  a  super¬ 
position  of  the  N^. 

CD 

M  «  ^  (49) 

n=0 


The  cylinder  func  tions  Z  p  of  a  noninteger  parameter  p  can  be  expressed 
in  given  polynomials  o.f  the  argument  and  sin  and  cos  functions.  For 
instance,  for  n  =  0, 


Z1  — y-  sin  j  k  r. 

2  TT  k  I- 


With  Pq  *  1,  and  sin  ^  k  r  ,)  sin  h  k  r  =■  i  |exp  (kr)  -  exp  (-kr)J  , 


Mo  »  Ao 

L  1-^ 

The  constant  factor  (2  i  V"?  j.s  hereby  combined  with  Aq* 


(50) 


Eeich  of  the  e-function  terms  in  the  bracket  of  Eq.  (50)  fulfills  the 
differential  equation.  Hence,  if  the  boundary  condition  requests  a  pole 
(Mq  CD )  at  the  origin  r  =  0,  only  the  second  term  is  used.  This  will  then 
lead  to  the  potentl.8l.  function  of  a  point  source. 


According  to  Eq,.  (27),  the  potential  function  is  given  by 

0  «  M  X"^, 
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Using  only  the  second  term  in  the  hracket  of  Eq,.  (50)#  and  letting  Aq 

take  care  of  the  sign  and  the  constant  factor  the  potential  function  of 

a  point  source  is  obtained. 

Mcp  [^k(r  t  z)] 


In  a  manner  very  similar  as  before,  it  is  possible  to  caiculate  from 
Eq,.  (51)  the  field  strength  in  the  r  direction  on  the  surface  of  a  very  small 
sphere  around  the  current  source  I,  convert  the  field  into  the  current 
density,  and  integrate  over  the  surface.  This  gives  the  strength  of  the  crir- 
rent  source  I,  and  the  constant  can  be  deteiminedl.  It  is 

whereby  denotes  the  conductivity  vsilue  at  the  current  source.  From 
Eqs.  (51)  eknd  (52)#  the  final  solution  is  obtained; 


For  k  *  0,  Eq.  (38)  is  X  k  exp  (2  ka)  ■  X^#  a  constant,  and 
Eq.  (53)  simplifies  to 


(‘-■rvTS-i- 

This  is  the  right  solution  for  the  point  source  in  a  medium  of  constant 
conductivity. 

The  equation  (53),  which  follows  here  as  the  Bla5)le8t  solution  of  Eqs. 
(48)  and  (49)#  is  the  key  solution  in  the  calculation  of  the  electric  field 
of  thunderstorms  of  Holzer  and  Saxon# ^  in  the  calculation  of  the  recovery 
curve  of  lightning  discharges  by  Tamura#^^  and  la  the  calculation  of  the 
thunderstorm  generator  by  the  author.^ 

It  is  mentioned  here  that  a  similar  procedure  to  that  used  for  polar 
coordinates  led  to  the  eigen  functions  of  the  differential  equation# 

div  grad  M  -  k^M  s  0  in  cylindrical  coordinates  z#  R#  0.  Again  assuming 
rotational  symmetry# 

^  exp  [^±(1  -  n2)2  ^  .  (54) 

Zq  is  the  cylinder  function  of  the  order  zero,  n  is  here  not  confined  to  an 
integer#  but  may  assume  any  value.  For  n  s  0, 

Mo  =  Ao  exp  (±kz). 
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1.  i  kz 

It  is  seen  that  the  root  of  the  conductivity  function  X2  =  Xq2  e  is  the 
eigen  function  of  the  order  zero  in  cylindrical  coordinates.  In  cartesian 
coordinates  the  eigen  functions  are  given  by 

exp  (±c^x)  exp  (te^y)  exp  (te^z).  (55) 

C3_j  C2j  and  are  arbitrary  functions  of  n,  restricted  only  by  the  equation 

(56) 

A  final  remark  is  made  about  the  representation  of  X  in  polar  coordi¬ 
nates.  It  is  seen  from  Eq.  (48)  and  from  Eq.  (50)  that  a  pure  e-function  is 
not  an  eigen  function  in  the  polar  coordinate  system.  Therefore  it  is  not 
possible  to  obtain  a  simple  solution  of  the  potential  function  0  for 
X  =  Xq  exp  (2  kr).  The  bjst  that  can  be  done  is  to  choose  the  foUoving 
function  for  X; 

exp  [2k  (r  -  r^,)] 

X  *  - - -  .  (57) 


In  this  case  ,. 

N  r.-  x^  •  Xq^  exp  (-k  Tq)  '  To  ’  ^  , 

which  is  an,  eigen  function  according  to  Eq.  (50).  Ao  is  given  by 
Aq  *  Xs  Tjj  exp  (-k  r^).  If  r^  denotes  the  earth’s  radius,  and  this  discus¬ 
sion  is  confined  to  the  space  between  the  earth’s  surface  and  about  100-km 
altitude,  then  tb.e  in  the  denominator  changes  very  little,  but  the 
e-functioa  in  the  r.oniinator  increases  in  the  requested  way.  Hence,  the 
representation  of  \  by  Eq.  (57)  absolutely  feasible.  This  will,  then, 
bring  agoin  the  full  benefit  of  the  sla^ile  solutions  as  outlined  above. 

The  Mirror  Lav  in  the  Current-Flow  Theory 

To  introduce  the  earth's  surface  as  an  equipotential  layer,  the  mirror 
law  is  applied  in  the  electrostatic  theory.  As  this  is  a  powerful  method 
which  leads  to  simple  solutions  of  boundary  problems,  its  application  in 
the  current- flow  theory  will  be  discussed  briefly. 


Take  a  point  source  I  and  place  it  in  a  cylindrical  coordinate  system 
z,  R,  ^  on  the  positive  7.  exis  at  the  point  z  »=  h.  With  regard  to  the 
coordinate  system  used  t'oi'  forniula  (53)^  fbe  zero  point  has  been  shifted 
down  the  z  axis  by  the  distance  h.  As  can  be  easily  verified,  this  means 
that  the  pcjtential  function  0  Is  now  expressed  in  the  following  form: 


exp  r  "k(r  +  z  ~  h)J 

^h  ""  \>  kh)  Is  here  the  conductivity  at  the  altitude  h,  and  r  is 

given  by  r  +  (-  "  h)'^  j 


(56) 
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To  now  introduce  the  plane  2  »  0  (eeirth  surface)  as  an  equipotential 
layer,  place,  analog  to  the  electrostatic  theory,  a  point  source  of  the 

strength  I*  at  the  mirror  point  of  I,  i.e.,  at  z  e  -h.  The  potential 
function  of  this  source  will  he 


I* 

f  TT  \* 


exp 


[^-k(r*  +  I  +  h)^ 


(59) 


with  r*  *  jjt  +  h)®  +  and  X*  B  Iq  exp  (-2k  h),  the  conductivity  value 
at  the  point  z  *  -h. 


The  superposition  of  and  resiilts  in  a  potential  function  0,  which 
will  he  zero  for  z  *  0. 


0  .  0j^  +  0* 


I 


exp  [-k(r  +  z  -  hj] 


^  **p  [-lt(r*  +  *  +  h)3 

4-  tr  X*  r* 

For  2  «  0  is  r  =  r*,  0  =  ~  exp  (kh)  +  ^  exp  (-kh),  or 


I*  »  - 

^h 


(60) 


(61) 


If  Eq.  (61)  is  inserted  in  Eq.  (60),  the  final  formula  for  the 
potential  function  0  is  obtained; 

*  '  [^a4lSEt.2S£j^]«p  -it-h)  .  (62) 

For  k  B  0,  Eq.  (62)  will  change  to  the  potential  fimction  0g  of  the 
electrostatic  theory. 

■  inn^  [i  -  3]  •  ‘^3) 

■With  X*  exp  (2  kh)  =  Xq,  from  Eq.  (61), 

I*  =  -  ^  I.  (64) 

Xh 

In  the  electrostatic  theory,  the  mirror  source  is  placed  at  the  mirror 
point  -h  and  its  strength  is  of  the  same  amotuit  hut  of  opposite  sign  of  the 
original  source.  It  is  seen  here  that  in  the  current-flow  theory  the  loca¬ 
tion  and  the  sign  reversal  are  retained,  but  that  the  amount  of  the  mirror 
source  is  smaller  than  the  original  source  in  the  same  proportion  as  the 
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conductivity  at  the  groxmd  is  smaller  compared  to  the  conductivity  Xji 
at  the  hei^t  of  the  origiaeil  source. 

Applied  to  the  current  flov  of  a  thunderstorm,  this  result  has  far- 
reaching  consequences,  as  pointed  out  by  the  author.®  It  means  that  only  a 
part  of  the  conduction  current  produced  by  the  thunderstorm — namely,  that 
given  by  I*— flows  to  the  earth’s  surface  in  the  Immediate  nei^borhood  of 
the  storm,  while  the  other  part  is  drained  off  to  the  ionosphere  and  con¬ 
tributes  to  the  air-eaxth  current  of  the  fair-weather  areas.  The  difference 
between  the  lines  of  current  flow  of  a  point  source  In  a  medium  of  constant 
conductivity,  X  «  (equivalent  to  the  field  lines  of  the  electrostatic 
case),  and  those  of  a  point  source  in  a  medium  of  increasing  conductivity, 

X  ■  exp  {2  h*),  can  be  easily  recognised  from  Figs.  3a  and  3b.®  This  dif¬ 
ference  also  affects  very  drastically  the  amount  of  charges  of  the  thunder¬ 
storm,  as  calculated  from  field  measurements  at  the  ground,^^^  ^5  the  field 
reversal  at  the  ground  from  a  bipolar  thunderstorm,  the  recovery  curve  of 
the  lightning  stroke, and  many  other  phenomena.  But,  as  the  object  of  this 
report  is  to  outline  the  methods  of  calculation  more  than  to  discuss  specific 
application,  the  mathematical  aspects  will  be  continued. 


The  Potential  Function  of  a  Deeaj 


Cfurrent  Source 


This  problem  finds  its  application  in  the  calculation  of  the  recovery 
curve  after  a  li^tnlng  discharge.  Field  records  of  a  thiuiderstorm  tadcen  at 
the  ground  show  that  each  li^tnlng  flash  increases  (or  decreases)  the  so- 
called  stationary  field  of  the  thunderstorm  very  suddenly,  and  that  after  the 
flash  is  completed  the  field  returns  to  its  preflash  value  in  an  e-function 
fashion,  whereby  the  time  constant  of  the  decay  Is  approximately  ten  seconds. 
The  retxim  is  called  the  recovery  curve. 


This  phenomenon  is  rather  puz2.ling  because  the  time  constant  of  the  air 
at  the  groimd  is  about  600  seconds,  and  charges  separated  or  produced  by  the 
lightning  flash  should  theref.ire  decay  nuich  slower.  Tamura^®  showed  in  an 
excellent  analysis  that  a  much  faster  decrease  of  the  field  at  the  ground 
will  result  if ,  besides  the  decay  of  the  lightning  charges  with  the  time 
constant  in  the  cloud  and  the  regenerating  effect  of  the  thunderstorm*  s 
charging  mechanism,  the  screening  effect  of  the  space  charge  is  taken  into 
account,  which  builds  up  under  the  influence  of  the  conductivity  which 
increases  with  altitude.  In  other  words,  his  calculations  are  based  on  the 
theory  of  current  flov.  The  calculation  given  here  may  be  considered  as  a 
part  of  Tamara's  analysis,  but  presented  from  a  different  point  of  view, 
namely,  to  show  a  classical  exaiiq)le  of  the  full  application  of  the  current- 
flow  theory. 


The  key  problem  is  that  of  the  decaying  point  charge.  Assume  that  the 
charge  is  deposited  by  the  lightning  flash  at  a  point  at  the  altitude  h. 
(if  it  is  a  cloud  flash,  in  addition  a  negative  charge  -Qq  wotild  have  to  be 
placed  at  a  higher  altitude  H  and  the  resulting  fields  superis^josed. )  The 
charge  decays  with  the  time  constant  9*  at  this  altitude,  gt  _ 

charge  Q  left  at  the  time  t  is  then  given  by  • 

-  A'i’. 

c 


Q  «  Qo  e 
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The  current  output  I  is  given  hy 

-  ^ 

I  =  -  ^  =  XL  e  ®  , 

dt  e  ^ 

Therefore,  the  tloie  function  T  of  the  c-urrent  source  is  given  hy  the 
escponential  function  of  Eq^.  (66). 

- 

e 

T  =  e 


The  conductivity  X  depends  here  only  on  the  space  coordinates,  hut  is  con¬ 
stant  with  time.  The  wei^ting  time  functions  Tg,  Tg,  and  T.  given  hy 
Eqs.  (17),  (18),  and  (19)  simplify  to 
_  Xt  \t 


Tg  =  e 


e  r  ^  . 

J  -rr  e  dt. 


- 

‘  I 


Xt 

e 

T  e  dt. 


To  »  e  « 


With  T  given  hy  Eq.  (67),  the  integrals  can  he  solved,  and  it  hecomes 


\  -  X’ 


Tg  = 

The  equation  Tg  +  Tg  =  T  is  fxilfilled  hy  Eqs.  (7I)  and  (72).  The  next  st^^p 
is  to  determine  from  the  initieO.  condition.  For  t  a  0,  it  Is  foiuid  that 

the  field  S’  shall  he  the  electrostatic  field  Fg,  since  the  lightning  flash 
occurs  in  such  a  short  time  that  the  space  charges  of  the  stationary  ctirrem  ' 
flow  field  due  to  the  conductivity  variation  have  not  accumulated.  From 
Eqs.  (21),  (67),  (70),  (71),  and  (72)  follows  for  1=0, 

^e  =  ^o  “  T'.rirr  Pg  +  x  -  X' 


"  X  -  X*  + 


— - —  7 
X  -  X*  ^s* 


If  Eq.  (73)  is  introduced  into  Eq.  (21),  obtained  with  regard  to  Eqs.  (71) 
and  (72)  is 
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(74) 


Eqiiatioi.  (74)  is  the  final  form  of  the  field  eqioation  of  the  decaying 
X)oint  source.  It  is  valid  for  every  point  in  space.  But  it  jnust  "be  reinein- 
benfrd  that  vector  addition  is  required,  because  Fg  and  Fg  do  not  have  the 

saane  direction.  Only  at  the  ground  is  the  direction  of  F,  and  the 
same,  namely,  vertical  to  the  ground  siirface.  In  this  case  the  vector 
equation  (74)  sijaplifieB  to  a  scalar  equation.  Notice  that  the  sum  of  the 

wei^ting  time  furxtlons  of  Fg  and  a^ain  is  T,  the  time  function  of  the 
current  source. 
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